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Abstract 

As a consequence of the wave nature of light, invisibility devices based 
on isotropic media cannot be perfect. The principal distortions of invisibility 
are due to reflections and time delays. Reflections can be made exponentially 
small for devices that are large in comparison with the wavelength of light. 
Time delays are unavoidable and will result in wave-front dislocations. This 
paper considers invisibility devices based on optical conformal mapping. The 
paper shows that the time delays do not depend on the directions and impact 
parameters of incident light rays, although the refractive-index profile of any 
conformal invisibility device is necessarily asymmetric. The distortions of im- 
ages are thus uniform, which reduces the risk of detection. The paper also 
shows how the ideas of invisibility devices are connected to the transmuta- 
tion of force, the stereographic projection and Escheresque tilings of the plane. 

PACS 42.15.-i, 02.40.Tt 
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1 Introduction 



The bending of light in dielectric media pQ is the cause of many optical illusions. 
For example, in a mirage in the desert |2J, light rays from the sky are bent above the 
hot sand where the air is thin and the refractive index is low. In this way the rays 
minimize their optical paths according to Fermat's Principle JT . They are creating 
images of the sky that deceive the observer as illusions of water [2] . Imagine a differ- 
ent situation [HHHE] where a medium guides light around a hole in it such that the 
rays leave the medium as if nothing were there, see Fig. ^ Any object placed inside 
would be hidden from sight. The medium would create the ultimate optical illusion: 
invisibility 2]. Recently, ideas for designing such invisibility devices have been dis- 
cussed jH Ej. Ideas for minuscule invisible bodies (smaller than the wavelength of 
light) are older but first schemes for implementations and interesting twists and 
methods have been developed recently H]. Cloaking devices require unusually 
strong refractive-index profiles, but it is conceivable that they can be built with di- 
electric metamaterials jZHHUHI- Such devices would operate in the microwave region 
of the electromagnetic spectrum [H] and perhaps also in some frequency windows in 
the visible range jl()j . 




Figure 1: Light propagation in a conformal invisibility device. The light rays are shown 
in red. The brightness of the green background indicates the refractive-index profile. The 
device is based on the conformal map ()A3|) and Maxwell's fish eye ()22l) with = 4ro and 
riQ = 2 (further details given in the Appendix). The device consists of an exterior and an 
interior layer with a clear boundary illustrated by the tiling in Fig. |3J The invisible region 
is shown in black. Anything could be placed there. The left figure illustrates how light is 
refracted at the boundary between the two layers and guided around the invisible region 
where it leaves the device as if nothing were there. In the right figure, light simply flows 
around the interior layer. 



2 



Strictly speaking, ideal invisibility devices based on isotropic media are impossi- 
ble due to the wave nature of light ^] . Highly anisotropic media, however, may 
lead, in principle, to the construction of perfect invisibility devices (4j. Expressed 
in mathematical terms, the inverse scattering problem for linear waves in isotropic 
media has unique solutions [TT]. Therefore, the asymptotic behavior of propagation 
through empty space, or a uniform medium, is only consistent with the actual prop- 
agation through a uniform medium. In theory, nothing can be hidden. In practice, 
a dielectric invisibility device would perhaps create a slight haze, instead of a per- 
fect image. The principal distortions of invisibility are due to reflections and time 
delays. Reflections can be made exponentially small for devices that are large in 
comparison with the wavelength of light jSj. Time delays are unavoidable. Unless 
the phase delays are tuned to be multiples of 2tt for some certain frequencies, they 
will result in wave-front dislocations at boundaries that lead to image distortions 
due to diffraction [Tj. Wave- front dislocations also pose the risk of detection by 
sensitive wave- front sensors jT3j. In this paper we calculate the time delay caused 
by the scheme jS] based on optical conformal mapping, assuming two-dimensional 
light propagation from infinity. We find that the delay is uniform for all directions 
and impact parameters, although the refractive- index profiles of invisibility devices 
are necessarily asymmetric ^3]. Therefore, the distortions of images composed of 
various spatial Fourier components are uniform, which reduces the risk of detection. 



2 Theory 

Our theory is based on geometrical optics pQ and in particular on Fermat's Principle 
and on Hamilton's analogy pQ between the propagation of light in media and the 
motion of particles in classical mechanics ^2]. Suppose that the refractive index 
profile n(r) does not vary much over scales comparable with the wavelength of light. 
In this regime of geometrical optics both polarization components of light ip for 
frequencies u obey the Helmholtz equation [T] 

V 2 + n 2 ^^ = (1) 

where c denotes the speed of light in vacuum. The Helmholtz equation (JTJ) is equiv- 
alent to the stationary Schrodinger equation with potential U and energy E such 
that [U 

U-E = -£. (2) 

Therefore we expect that Hamilton's equations for light rays are equivalent to New- 
ton's equations of mechanical particles 1 moving in the potential (J2J) • The frequency u 
plays the role of the Hamiltonian and the wavevector k corresponds to the canonical 
momentum, 

ck 

uj = — , k = Ikl . (3) 
n 

1 Light rays in moving media behave like particles in magnetic fields for low velocities |16U17llT§| 

and like particles in gravitational fields ^3 ^1 120] in general. 
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Indeed, we obtain from Hamilton's equations ^2] the relations 

dr dui c k c , dk dui ck _ , c<j _ 9 . „ . 

17 = ^- = -T = ^k, — = - — = Vn 2 = — Vn 2 (4) 
at ok n k n z u dt or ln A 2n z 

that result in the equation of motion for light rays 

n 2 d n 2 dr Vn 2 



c dt c dt 2 
We can express this equation as Newton's second law 

d 2 r Vn 2 



(5) 



(6) 



dr 2 2 

with the effective time increment dr measured in spatial units and defined by 

cdt = n 2 dr. (7) 

Equation (jHJ) also reveals the connection to Fermat's Principle pQ: light in media 
with refractive index n takes the shortest (or longest) optical path where the optical 
path length is defined, in Cartesian coordinates, as 



/ tia/ da; 2 + dy 2 + dz 2 



To see this we use the fact that the modulus of the Hamiltonian velocity v equals 
c/n, a simple consequence of Hamilton's equations and write 

n d n dr Vn 2 

v dt v dt ~ ~2~ ' ^ ' 

These are the Euler-Lagrange equations ^3] for the effective Lagrangian nv. Hence 
they minimize or maximize the integral of nv that is equal to the optical path length 
(JHJ), which proves Fermat's Principle. The phase of a light ray is given by |T] 

k-dr-uot. (10) 

Along a ray trajectory the phase is constant. Consequently, the phase delay J k-dr 
corresponds to uot. Therefore, the Hamiltonian time t measures the true time delay 
of light caused by the refractive- index profile, whereas the Newtonian time r serves 
as a convenient parameter to characterize the ray trajectories. 

Another ingredient of our theory is optical conformal mapping [Sj. Consider 
an effectively two-dimensional case where the medium is uniform in one direction 
and the light propagates in a plane orthogonal to this axis. It is convenient to use 
complex numbers z = x + iy for describing the Cartesian coordinates x and y in this 
plane. In complex notation, the Helmholtz equation assumes the form 

4 Jk + 4)^°- < u > 



4 



Figure 2: Optical conformal map. A dielectric medium conformally maps physical space 
described by the points z = x + iy of the complex plane onto a stack of Riemann sheets 
if the refractive-index profile is \dw/dz\ with some analytic function w(z). An invisibility 
device [3] consists of two layers in real space, as Fig. ^ indicates. On the Riemann surface, 
the top sheet corresponds to the exterior and the first lower sheet to the interior layer. The 
figure illustrates the typical fates of light rays in such media. On the Riemann sheets rays 
propagate along straight lines. The rays shown in blue and green avoid the branch cut 
and hence the interior of the device. The ray shown in red crosses the cut and passes onto 
the first lower sheet where it approaches oo. However, this oo corresponds to a singularity 
of the refractive index and not to the oo of physical space. Rays like this one would be 
absorbed, unless they are guided back to the exterior sheet. 



Suppose that the complex z coordinates are transformed to new coordinates w with 
an analytic function w(z) that does not depend on the complex conjugate z* and 
hence satisfies the Cauchy- Riemann differential equations [21]. Such analytic func- 
tions define conformal mappings of the complex plane onto Riemann surfaces [22], 
see, for example, Fig. [2] Since d/dz = (dw/dz) d/dw and d/dz* = (dw*/dz*) d/dw* 
we obtain in w space the Helmholtz equation (fTTj) with the transformed refractive- 
index profile n' that is related to n as I2H] 



n = n' 



dw 



dz 



(12) 



On the Riemann surface, light thus propagates according to the refractive index 
profile n' . The strategy [3] for designing an invisibility device is to take advantage 
of the sheets of the Riemann surface. Each sheet corresponds to a distinct region 
in physical space, see, for example, Fig. El The branch cuts of the Riemann surface 
represent the boundaries between the various regions. If one wishes to hide an object, 
one should hide it on Riemann sheets and prevent light from entering these sheets. 
To do this, light that has ventured across a branch cut into the interior of the device 
should be guided back to the exterior. This is done by placing a refractive-index 
profile on the first interior sheet in which all ray trajectories are closed [S]. The 
device thus consists of two layers, an outer layer that corresponds to the exterior 
sheet on the Riemann surface in w space and an inner layer that corresponds to 
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the first interior sheet. For the outer layer we require that the refractive index 
approaches unity at oo, the value for empty space, which implies for the conformal 
map 



At the boundary between the exterior and the first interior sheet light is refracted 
PP according to Snell's law 2 , unless it is totally reflected [TJ|3]. Since refraction is 
reversible, the light rays are refracted back to the original direction of incidence 
when they leave the branch cut to the exterior sheet. Seen on the Riemann surface, 
light rays perform loops around a branch point that guide them back to the exterior 
sheet, the outside layer of the device. Seen in physical space, light is guided around 
the object and leaves the cloaking layers of the device as if nothing were there. 



Figure 3: Riemann sheets are tiles. In optical conformal mapping, Riemann sheets 
represent regions of two-dimensional space, tiles of various forms |21l \22\ . The figure 
illustrates the tiling behind the light propagation of Fig. ^ (details given in the Appendix). 
The upper imaginary half-plane of each sheet corresponds to a grey tile and the lower half- 
plane to a white tile. The exterior and the interior sheets of the invisibility device occupy 
two of such pairs of tiles and the hidden core takes the rest. 



2 The law of refraction was discovered by the Arabian scientist Ibn Sahl more than a millennium 
ago PI]. 




for z 



oo . 



(13) 
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3 Time delay 



The time delay caused by an invisibility device depends of course on its spatial 
extension. Imagine that the refractive- index profile n(r) is replaced by n(£r) with 
the constant scale £. If r(i) and k(t) are solutions of Hamilton's ray equations (j4*|) 
then £r and k are solutions, too, if t is replaced by £t. This proves that the time 
delay is directly proportional to the spatial extension of the refractive-index profile, 
as one would expect. 

Since the conformal mapping w(z) is simply a coordinate transformation, the 
time delay between two points in z space and the delay between the corresponding 
points in w space are identical. Therefore, the part of the refractive-index profile n 
in physical space that is due to the optical conformal mapping does not influence 
the time delay at all. Delays are only caused by the index profile in w space that 
serves to guide light around a branch cut on the first interior sheet. For simplifying 
the notation, we denote this profile by n (dropping the prime). 

Suppose that the transformed refractive-index profile n on the Riemann sheet 
is radially symmetric with respect to one branch point and designed such that all 
trajectories of light rays are closed curves around that point. The time delay to of 
the invisibility device is equal to the time light takes to perform a loop around the 
branch point. A branch point where v sheets meet requires v turns iy is the winding 
number). To calculate the delay, we use polar coordinates r and ip centered at the 
branch point. We obtain from the conservation law of energy [15J for the Newtonian 
dynamics (JBJ) of light rays 

(^) 2 + r 2 ( ( ¥) 2 = 2(E-U)=n\r). (14) 



\dr J \dr 

As in the standard theory of motion in central potentials we also use the con- 
servation of the angular momentum, 

written here in terms of the impact parameter b, and obtain from the conservation 
laws jUD and jEJ) the relation [23] 

dp> = / dr . (16) 
ryn 2 r 2 — b 2 

To calculate the time delay, we express the Hamiltonian time increment ((ZJ) in terms 
of dcp, utilizing the conservation of the angular momentum ()15j) . Then we use the 
relation (fTo^) to write the time delay as an integral over the radial range of the 
trajectory. The range of r is bounded by the radial turning points r± where dr/dip 
vanishes, which implies 

n 2 (r ± )r| = 6 2 . (17) 

One trajectory between r_ and r + corresponds to half a turn around the branch 
point. Consequently, 

^= P / rdr = fV^^ + ^. (18) 
2u J r ^Jn 2 r 2 - b 2 J r _ r * v ' 
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In general, the time delay depends on the impact parameter. However, for closed 
loops, to turns out to be independent of b. To see this, we differentiate the time 
delay (|18j) with respect to the impact parameter b and obtain from Eqs. and 

(ED 

because, when all loops around the branch point are closed, reaches 7r regardless 
of the value of b. Consequently, the time delay does not depend on the impact 
parameter at which light has entered the branch cut to the first interior sheet, i.e. 
the interior layer of the device. The invisibility device causes a uniform time delay. 

4 Examples 

Reference |S] mentions two examples of refractive- index profiles on the interior sheet 
that can be used to circumnavigate the branch point such that all loops are closed, 
the harmonic-oscillator profile 



^1 -r 2 /rf (20) 



that is related to a Luneburg lens [2H1 EH] and the Kepler profile 



ri2 



y/r/r 2 ~ 1 (21) 



that is related to an Eaton lens [213 EH1 EZj- Here r% and r 2 are constants that 
describe the boundaries of the refractive-index profiles where n vanishes. Seen on 
the Riemann surface, light cannot penetrate the outside of circles of radii r% and 
r 2 , respectively, because here the refractive index would be purely imaginary. The 
optical conformal mapping turns areas on the first interior sheet inside out |2*T) . 
Therefore, the exterior of these circles corresponds to the invisible interior of the 
device. 

The harmonic-oscillator and the Kepler potential are the only spherically sym- 
metric potentials U where the trajectories for all bound-state energies E are closed 
|15j . However, what matters in the propagation of light rays is the difference (J2J) 
between U and E. Therefore it is sufficient when for a specific value of E the tra- 
jectories for all angular-momenta b are closed. A known example where this is the 
case is Maxwell's fish eye 123 E3 EE] with the refractive-index profile 

" , = i+W' (22) 

The constant radius r 3 characterizes the scale of the index profile, the radial half- 
point, and no is a parameter that defines the refractive index at the branch point. If 
Maxwell's fish eye is employed to guide light back to the exterior Riemann sheet, the 
entire interior sheet is reached by the incident light. However, when the Riemann 
surface contains more sheets than the exterior and the first interior sheet, all the 
remaining sheets are hidden. Anything placed there is invisible. 
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The Newtonian equation of motion for rays generated by the harmonic-oscillator 
profile (|2(Jj) describes Hooke's law of a force proportional to the distance, whereas 
the Kepler profile generates Newton's inverse-square law. In both cases, the 
trajectories form ellipses for all bound states, a fact that Newton found exceptionally 
remarkable [23 EDI- However, Hooke's law and Newton's law can be transformed into 
each other 3 by a transmutation of force according to the Arnol'd-Kasner theorem 

[3UE2E3]. 



4.1 Hooke's force 

In the case of the harmonic-oscillator profile (|20|) the ray trajectories are very simple 
[T*"] : they form a set of ellipses centered at the origin (the branch point), see Fig. 
["{J To calculate the time spent along a given ellipse we use Cartesian coordinates x 
and y rotated such that they match the axes of this ellipse. We describe the ellipse 
as 

x = acos£, y = 6sin£, C = T / r i (23) 

with the constants a and b being the axis lengths. One easily verifies that the 
trajectory ("2~3~|) solves the Newtonian equation of motion (0). The ray trajectory 
corresponds to the Newtonian energy 1/2, which implies 



Consequently, we obtain for the time delay 

"2tt / ^2 



r r 2 ^ f a 2 b 2 \ 
ct = j n 2 dr = vr\ J ( 1 j cos2 £ 2 s * n2 £ ) ^£ = z/7rr i • (25) 

In agreement with our general results, the time delay caused by the harmonic- 
oscillator profile depends on the spatial extension r x and is otherwise uniform. 



4.2 Newton's force 

In order to calculate the time delay caused by the Kepler profile (l2*Tj) . see Fig. we 
use the transmutation of force p?T) 13*2*1 13*3*] to the harmonic-oscillator profile, i.e. the 
transformation of Newton's law into Hooke's law that is also based on conformal 
mapping [22] ■ Consider trajectories in the complex plane, say the z plane, although 
in our case this plane is one of the Riemann sheets generated by the optical conformal 
mapping in the first place. Suppose that the trajectories are conformally mapped 
by the analytic function w(z). We obtain from Eqs. (J2J) and (fT!2|) the relations 



U-E = - 



7r 



n 



/2 



dw 



dz 



([/' - E') 



dw 



dz 



(26) 



3 Ironically, despite Newton and Hooke reportedly having been bitter rivals, their most celebrated 
force laws are essential identical ["""J [33 HS1 ■ 
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Figure 4: Light guiding using Hooke's law (similar to a Luneburg lens) with the 
refractive-index profile (|2U|) in w space. The device guides light that has entered its inte- 
rior layer back to the exterior, represented in the left figure using two Riemann sheets that 
correspond to the two layers, seen from above. The right figure shows the corresponding 
ray propagation in physical space with the optical conformal map l)A3|) and n = 8ro. At 
the branch cut in the left figure, the thick black line between the two points, the branch 
points, light passes from the exterior to the interior sheet. Here light is refracted accord- 
ing to Snell's law, see Eq. (jA8|) of the Appendix. On the lower sheet, the refractive-index 
profile (|2Uj) guides the rays to the exterior sheet in elliptic orbits with one branch point 
in the centre. Finally, the rays are refracted back to their original directions and leave on 
the exterior sheet as if nothing had happened. The dotted circle in the figure indicates 
the maximal elongations of the ellipses. This circle limits the region in the interior of the 
device that light does not enter. The outside of the circle and the other Riemann sheets of 
the map (|A<3|) correspond to the inside of the device in physical space, as the right figure 
shows. Anything inside this area is invisible. 



Consequently, if the potential U can be written as the modulus square of an analytic 
function, the potential U' is proportional to the modulus square of the inverse of 
this function expressed in terms of the new coordinates, 



U{z) = -E' 



dw(z) 



dz 



U\w) 



-E 



dz{w) 



dw 



w 



2n 



Ef=- 1 -. 
2 



(27) 



The trajectories are mapped onto each other by the transformation w(z). The 
potentials U and U' are thus related to each other, generating dual forces p?T] . 
Consider 



:28i 



The map w(z) corresponds to the Hooke potential (j2j) of the harmonic-oscillator 
profile (f2T)|) with energy E = 1/2 and, in turn, w(z) generates the Kepler profile 
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Figure 5: Light guiding using Newton's law (related to an Eaton lens) similar to Fig. 
131 with the refractive-index profile (|21|) in w space and the optical conformal map (jA3|) 
where r 2 = 8ro- The left figure shows how light is guided on a Riemann sheet, whereas 
the right figures shows the ray propagation in physical space. On the Riemann sheet, light 
rays form elliptic orbits with one branch point in the focal point, instead of the centre, as 
in Hooke's case. 



(|2*Tf that corresponds to Newton's inverse square law with the parameter 



Since conformal mapping does not influence the time delay in light propagation, 
the delay generated by the Kepler profile ()21j) corresponds to that of the harmonic- 
oscillator profile (|20jl . apart from one subtlety: already a half-ellipse of the harmonic- 
oscillator is mapped onto a complete Kepler ellipse, because of the square map ()28jl . 
Consequently, the time delay is 



in complete analogy to the result (J25|) for the harmonic-oscillator profile. The time 
delays are thus identical for identical ranges T\ and r 2 of the refractive-index profiles. 

4.3 Maxwell's fish eye 

Maxwell's fish eye turns out to represent another classic conformal mapping [23] > the 
stereographic projection discovered by Ptolemy and applied in the Mercator map 
projection 4 . Figure |H1 shows how the points on the surface of a sphere (X,Y, Z) are 
mapped onto a plane, say the z plane, according to the formulas [2U EH] 




(29) 



ct = -ivri = v-nr 2 , 



(30) 



1 • X -\- lY 2 . \f2 1 rv2 2 

z = x + ty= _ , X + Y + Z = r 3 



(31) 



4 The Mercator map projection is the logarithm of the stereographic projection 
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Figure 6: Stereographic projection, mapping the (x, y) plane onto the (X, Y, Z) surface 
of a sphere. A line drawn from the North Pole of the sphere to (x, y) cuts the surface of 
the sphere at (X,Y,Z). Circles on the plane are mapped into circles on the sphere and 
vice versa [35] . 



with the inverse 

x+iy = r$w> z = ^wtt- rW + ^ (32) 

We find that the square of the optical-length element (jSJ) for Maxwell's fish eye (122)) 
is 

ds2 = "» jtwS? = n i (dA ' 2 + dy2 + dz2) • (33) 

Consequently, the light rays of the fish eye are mapped into rays on the surface of 
a sphere with radius r 3 and uniform refractive index n /2. According to Fermat's 
Principle PQ, the rays are geodesies, lines of extremal optical path length, see Fig. 
[7| On a uniform sphere, the geodesies are the great circles. Since they are closed 
curves on the sphere, the light rays are closed on the plane as well, as we required. 
Furthermore, the stereographic projection maps circles onto circles [33] and so the 
light rays in Maxwell's fish eye ()22|) form circles in the plane PQ 1231 12H] ■ 

The calculation of the time delay of light circling in Maxwell's fish eye is el- 
ementary now, because to does not depend on conformal transformations and in 
particular on the stereographic projection (|3*T|): to simply is the time delay of light 
during v loops on the surface of a sphere with radius r 3 and refractive index no/2, 
which gives 

cto = v y 2vrr 3 = z/vrr 3 n . (34) 

In agreement with our general results, the time delay is uniform and proportional 
to the length scale of the refractive-index profile. 
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Figure 7: Light guiding using Maxwell's fish eye. The interior layer of the invisibility 
device is represented by a sphere of radius r% using the inverse stereographic map Q32JI . 
At the boundary, the branch cut, the light drops onto the sphere where it propagates in 
great circles. After jumping up to the exterior sheet the light rays leave the device. This 
behavior is generated by the refractive-index profile of Maxwell's fish eye (1221) that repre- 
sents the stereographic projection illustrated in Fig. El as an optical conformal mapping. 
The pictures illustrates light propagation for angles of incidence of ±7r/7. 

5 Conclusions 

In isotropic media, no illusion is perfect due to the wave nature of light jTTj. Con- 
sequently, conformal invisibility devices [5] cannot be perfect; they cause reflections 
and time delays. However, the reflectivity can be made exponentially small for 
macroscopic devices and the time delay is uniform for all directions and impact 
parameters. This is important, because images consist of light propagating in a 
range of directions, having a range of spatial Fourier components. The time delay 
occurs when the light reaches the interior layer of the device. It will cause wavefront 
dislocations at the two sides, unless the phase delay is tuned to be a multiple of 2n 
for some certain frequencies. The diffraction of light will slightly blur the image pQ , 
but the haze caused is uniform. 
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Appendix 



Most pictures of this paper are based on a conformal map using elliptic modular 
functions [33]. These functions are connected to many branches of mathematics, 
including the proof of Fermat's Last Theorem [33]. We adopt the notation of the 
Bateman Manuscript Project on Higher Transcendental Functions [B3] (not Neharis 
notation 5 ). We use the modular function J known as the Klein invariant, illustrated 
in Fig. [H] and expressed here in terms of the modular function A as 



4 (1-A + A 2 ) 3 
J( " ) = 27 A*(l-A)* ' X = 16q 



I oo \ 

^ g m ( m +!) 

m=0 



q" 



(Al) 



m=l 



J 



with 



q = e*" z , \q\<l, (A2) 
see Eqs. (36-38) of Ref. [S3]- Note that the expression for A is rapidly converging 



2 




-2-10 1 2 



Figure 8: The Klein invariant J{z) tiles the upper half z plane with an infinite sequence 
of circular arches. The arches continue near the real axis in an infinitely intricate structure 
(not shown here) . Identifying all horizontal strips of length 1 and deforming the outer arch 
to a circle leads to Fig. [3] that illustrates the tiling of the optical conformal map (|A3|) . 



and therefore well suited for numerical computations. Consider the map 
/ ln(432z/r ) ^ 31r 16r (1 - A + A 2 ) 3 31r 

w = 4roJ i 2^)-^ = ^T "IT (A3) 



with 



q = . . (A4) 

V4322/ro 



5 In Nehari's book (221 ^{ z ) is denoted as J{z). 
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The constant r characterizes the spatial scale of the optical conformal mapping. Far 
away from the device light should propagate through empty space, which implies 
that w ~ z for z — > oo. The map (jA3|) is chosen such that this is the case. Indeed, 
we obtain from the representation (jAl|) the first terms of the Laurent expansion 



1 31 1823 2 . 

J ~ 9 + ^ + ? fOT 9^°) ( A5 ) 



1728g 2 72 16 
which implies, according to Eqs. (|A3|) and (jA4 



!823 r 2 , r 

~ z + — for z — > oo . (A6) 
1728 z 



In the exterior of the device, the map (|A3|) approaches the simple example considered 
in Ref. whereas in the interior the map (|A3|) represents the infinitely more 
complicated Riemann surface illustrated in Figs. El and |H1 The Riemann surface 
contains three branch points [22 El], wi = (41/18)ro, w 2 = — (31/18)r and = 
00 with winding numbers 1, 2 and 00, apart from the exterior sheet with only Wi 
and w 2 [21 Oil- 
In order to calculate the ray trajectories, we consider the ray dynamics in w 
space and then transform w to the physical trajectories in z space. We describe 
both the ray trajectories w and the wavevectors k by complex numbers. In the 
exterior sheet light propagates along straight lines. Given a point w on the exterior 
sheet, we numerically solve Eq. (jA3|) for z using the inversion z(w) of the asymptotic 
map (|A6|) as the starting value, 

z ~ I2 ( 36w + v 7 ^ 432 ^ 2 ~ 1823 )) ( AT ) 

for Imu; > and we use z*(w*) otherwise. At the branch cut between the exterior 
and the first interior sheet light is refracted according to Snell's law PjJ 121] . Since 
the modulus of k equals n(u/c), we obtain for a light ray incident at the angle tp 
Snell's law in complex notation as 



c 



k = — ( sin if — iyjn 2 — sin 2 ip ) . ( A8) 



On the first interior sheet, we solve Hamilton's equations for a radially symmetric 
index profile around the branch point W\, 

dw k dk n r (\w — Wi\) \k\ w dn(r) 

dl n(\w — wi\) \k\ ' dZ n 2 (\w — w x \) \w\ ' dr 



for the propagation distance I = ct and Maxwell's fish eye (|22j) with n = n 3 , r 3 = 4r 
and no = 2. The parameters are designed such that the refractive index on the 
Riemann structure reaches 1 at the other branch point w 2 and exceeds 1 along the 
branch cut. In this way, total reflection is excluded for the lowest possible value of 
Uq. To calculate z(w) on the first interior sheet we utilize the modular symmetry of 
the Klein invariant |34j . 

J(z') = J(z) for z = -z- 1 , (A10) 
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which leads to 

^ = — e -p ( , Z" 2 / ^ ( An ) 

432 P V ln ( 432 ^/ r o) / 

for the position 2' that shares the same numerical value of w as z, but corresponds 
to the first interior sheet. We follow the same procedure as for the exterior layer of 
the device to calculate z and then transform it into z' to continue the trajectory in 
the interior layer. 
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